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l"~5 ■ Abstract. We study the evolution of a partonic jet shower propagating through a 

quark-gluon plasma. Combining the in-medium evolutions of the leading parton and 

shower gluons, we compute the depletion of the energy from the jet cone by dissipation 

^3 • through elastic collisions with medium constituents, by scattering of shower partons to 

Mh| larger angles, and by radiation outside the jet cone. Numerical results are presented 

O i' for the nuclear modification of dijet asymmetry in Pb+Pb collisions at the LHC. 

1. Introduction 

t— i ; 

Jet quenching has been studied in some details (e.g., see [Ul2]) in relativisitic heavy-ion 
collisions due to the fact that hard partons interact with the traversed mater and lose 
some of the initial energy by elastic and inelastic collisions with medium constituents, 
providing powerful probes to the hot and dense matter produced in these reactions 
[31 IH 0, E] . One of the most spectacular jet measurements from the first heavy- ion runs 
at the CERN large Hardon Collider (LHC) is the nuclear modification of the energy 
^ asymmetry between two correlated jets emitted in opposite azimuthal direction around 

c<3 ■ the beam axis in Pb+Pb collisions at y/s^N = 2.76 TeV [7J |8]. 

In this work, we study the medium modification of a parton shower propagating 
through a quark-gluon plasma [9]. As the jet shower traverses the medium, we 
calculate the energy loss experienced by the jet shower within a cone angle defined 



as R = J(A(j)) 2 + (A77) 2 , which includes collisional energy loss by the leading parton 
and shower gluons through elastic collisions with medium constituents, as well as the 
energies carried by shower gluons that are radiated or scattered outside the jet cone. 
We present the numerical results for the nuclear modification of dijet energy asymmetry 
in Pb+Pb collisions at the LHC. 

2. Jet shower evolution in medium 

In our model, a jet shower is described by two quantities: the energy El of the leading 
parton and the double-differential distribution f g (u,k±,t) = dN g (u,k±,t)/du)dk\ of 
shower gluons, where u denotes the gluon energy and k± the transverse momentum 
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with respect to the jet axis. For a jet defined by a cone angle R, the energy contained 
in the jet cone is given by 

Ej(R) =E L + E g (R) = E L + I iodudklf g (oo, k\), (1) 

J Ft 

where the subscript R denotes the integration taken within the jet cone, k± < uR. 

The evolution of a jet shower propagating through medium is governed by two 
equations, i.e., the evolutions of the leading parton energy and shower gluon distribution: 

r r dN med 

E L (t f ) = E L (ti) - J e L (t)dt - J ududk 2 ± dt dud 9 k 2 dt > ( 2 ) 

d t( * *\ ~ d f9 ^Y72 f dN T CA /ox 

Here we include both the contributions from elastic collisions and medium- induced 
radiation, as well as that from transverse momentum broadening experienced by shower 
gluons. Note that a source term from vacuum radiation should be also included for the 
interference with medium-induced radiation. In this application, we assume that jets 
experiences vacuum radiation before interacting with medium, thus vacuum radiation 
serves as the initial conditions for jet shower in-medium evolution. 

After solving the above jet shower evolution equations, one may calculate the total 
energy loss from the jet cone, 

AEj = Ej(ti,R)-Ej(tf,R). (4) 

A few inputs must be supplemented to solve the above jet shower evolution 
equations and calculate the energy loss from the jet cone. The initial gluon distribution, 
f g (u,k\,ti), before interacting with the thermal medium at time U, is generated from 
PYTHIA [10]. We impose a formation time cut-off, i.e., only gluons with a formation 
time Tf = 2Ex(l — x)/k\ smaller than U are radiated, where x = u/E denotes the 
energy fraction of the radiated gluons. The rate for the medium-induced gluon radiation 
is taken from the higher-twist formalism of jet quenching [TTJ | 

dN™ ed 2a s xP(x)q{t) . 2 t-t 



dudk\dt 7r uk\_ 2tj 

where P(x) is the vacuum splitting function. We further relate two transport coefficients 
by the fluctuation-dissipation theorem, q = 4Te, assuming that the medium is in thermal 
equilibrium. When solving for the radiated gluon distribution f g (oj,kj_,t), a lower cut-off 
of 2 GeV is imposed on the energy of shower gluons. 

3. Nuclear modification of dijet asymmetry at the LHC 

We now calculate the medium modification of the dijet asymmetry in Pb+Pb collisions 
at the LHC. The space-time profile of the medium temperature T(f±,t) at midrapidity 
is modelled as follows. The initial entropy density s ~ T 3 of the medium is set 
to be proportional to the density of participating nucleons in the colliding nuclei, 
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with the nuclear density distributions taken as Woods-Saxon profiles. The medium 
created in Pb+Pb collisions at ^snn = 2.76 TeV is assumed to thermalize at 
to = 0.6 fm/c, at which the temperature of the hottest point in central collisions 
is set to be T = 520 MeV. The time evolution of medium is modeled by a one- 
dimensional boost-invariant expansion, i.e., the temperature falls with time as i" 1 / 3 . 
For jet shower evolution, we assume that it experiences vacuum radiation before to, after 
which jet-medium interaction starts until the local temperature of the medium drops 
below T c = 160 MeV. The transport coefficients are scaled according to the temperature 
of the medium, q ex T 3 , with a constant factor adjusted to fit the experimental data. 
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Figure 1. (Color online) Distribution of dijet energy asymmetry factor Aj for p + p 
and Pb+Pb collisions at ^/sjvjv = 2.76 TeV at the LHC. Left two panels: 0-10% and 
10-20% centralities compared to the ATLAS [7[. Right two panels: 0-10% and 10-20% 
centralities compared to CMS [8]. 



The results of the medium-modification of dijet asymmetry factor Aj is shown in 

Fig. CEJ where Aj is defined as 

E T ,i - E T2 , „v 

Aj = F I F ' (6) 

with E T>i , (i = 1, 2) the transverse energy of the leading and sub- leading jet, respectively. 
The ATLAS Collaboration measured this quantity with the trigger jet Ex,i > 100 GeV 
and the second jet in the opposite hemisphere A0 > 7r/2 with Et,2 > 25 GeV. For CMS 
Collaboration, these values are Et,i > 120 GeV, A(f> > 2n /3, and Et^ > 25 GeV. 

The vacuum dijet events for p+p collisions at the LHC energies shown by the shaded 
yellow mountains are generated from PYTHIA [10], where jets are reconstructed using 
anti-^T algorithms [IB]. Note a Gaussian smearing with width oc \/E~~j is applied here to 
taken into account the detector response and other smearings; this reduces some amount 
of very symmetric dijet events. The modification of each dijet event in Pb+Pb collisions 
is obtained as follows. For each dijet event, its production point is sampled according to 
the distribution of the binary nucleon-nucleon collisions from colliding two lead nuclei. 
Then additional energy loss from the jet cone due to jet-medium interaction is applied 
and the dijet asymmetry Aj distribution is recalculated from the surviving dijet events 
as shown by the thick green lines. The trigger bias is approximated by letting the 
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leading jet propagate along the shorter path, and the subleading jet propagate along 
the other direction when dijets are asymmetric (Aj > 0.1). For nearly symmetric jet 
pairs (A j < 0.1), we do not apply such trigger bias treatment. 

As expected, the energy asymmetry of dijes is significantly increased by the in- 
medium evolution and more prominent in the most central Pb+Pb collisions than in less 
central events. By fitting to the data we obtain the length-averaged transport coefficient 
in central collisions (q) ~ 0.9 GeV 2 /fm (with about 20% difference between ATLAS and 
CMS for the best descriptions of the data). This corresponds to q ~ 2 GeV 2 /fm at 
T = 400 MeV, the highest temperature available in Au+Au collisions at Relativistic 
Heavy Ion Collider (RHIC). 

4. Summary 

In summary, we have presented a simplified but realistic model for studying the evolution 
of a jet shower propagating in a quark-gluon plasma. We have applied the model 
to compute the nuclear modification of dijet energy asymmetry in Pb+Pb collisions 
at the LHC. The observed dijet asymmetry can be described with values of parton 
transport coefficients similar to those describing jet quenching data at RHIC. Some 
further directions include a complete Monte-Carlo simulation of jet shower evolution in 
medium and reconstruction, and realistic event-by-event hydrodynamical simulation of 
fluctuating background (e.g., see [H]) and its effect on dijet energy asymmetry [1~5] . 
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